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Abstract 

On the proper timescale, amorphous solids can flow. Solid flow can be ob- 
served macroscopically in glaciers or lead pipes, but it can also be artificially 
enhanced by creating defects. Ion Beam Sputtering (IBS) is a technique in 
which ions with energies in the 0.1 to 10 keV range impact against a solid tar- 
get inducing defect creation and dynamics, and eroding its surface leading 
to formation of ordered nanostructures. Despite its technological interest, 
a basic understanding of nanopattern formation processes occurring under 
IBS of amorphizable targets has not been clearly established, recent exper- 
iments on Si having largely questioned knowledge accumulated during the 
last two decades. A number of interfacial equations have been proposed in 
the past to describe these phenomena, typically by adding together different 
contributions coming from surface diffusion, ion sputtering or mass redistri- 
bution, etc. in a non-systematic way. Here, we exploit the general idea of 
solids flowing due to ion impacts in order to establish a general framework 
into which different mechanisms (such as viscous flow, stress, diffusion, or 
sputtering) can be incorporated, under generic physical conservation laws. 
As opposed to formulating phenomenological interfacial equations, this ap- 
proach allows to assess systematically the relevance and interplay of different 
physical mechanisms influencing surface pattern formation by IBS. 

Keywords: Ion Beam Sputtering, Hydrodynamics, Solid flow, Pattern 
formation, Stability, Viscous flow, Stress, Erosion 
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1. Introduction 



Observations of nano-scale patterns on the surfaces of solid targets that 
undergo ion-beam sputtering (IBS) by 100 eV to 10 keV ions, date back at 
least to the early 1960's [H, HJ. They correspond either to amorphous mate- 
rials like glass, or to targets that are amorphized by this type of irradiation, 
like semiconductors [i[ • The fascinating resemblance with macroscopic struc- 
tures, like ripples on water or sandy dunes0 emphasized in j^], goes beyond 
a mere visual analogy 0, IE[ , underscoring the basic interest of this type of 
structures. Moreover, the high efficiency of IBS to produce surface nanopat- 
terns (ripples and dots) over large areas on top of a wide variety of targets 
including metals and insulators p, [(J , furnishes the technique with a large 
potential for applications. 

Specifically, in IBS both the erosive process and the appearance of the 
patterns are due to the complex ion-target interactions: thus, ions create per- 
manent defects in a surface layer whose thickness is of the order of the aver- 
age penetration depth, trigger material ejection (sputtering) from its surface 
and induce material redistribution and flow jsl, [if. A conceptual framework 
explaining pattern formation in amorphous or amorphizable systems is avail- 
able since the late 80's, through the work by Bradley and Harper (BH) [l(| 
that was based on Sigmund's theory of linear collision cascades. Physically, 



a morphological instability [ll| originates in the higher erosion rate (yield) 
at surface minima as compared with surface maxima, in such a way that a 
pattern forms for any value of the angle of incidence and ion energy. This 
has become the theoretical paradigm supporting most experimental observa- 
tions. Minor discrepancies with this view have been attributed to secondary 
effects, leading to refinements of BH's theory (see reviews in jEl, 0]). 

However, many of the experiments done on elemental targets have been 
shown to be affected by contamination [jjj], suggesting a correlation between 
pattern formation and a sizeable concentration of undesired species. In view 
of this, very recently clean experiments have been specifically tailored in order 



to avoid this feature |!3Hl6|. taking Si targets as representative cases for the 



large class of substrates that become amorphous under irradiation. Remark- 



In Fig. Q] we show the nice naked-eye similarity between water and Silicon ripples at 
very different orders of magnitude. 
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Figure 1: Water ripples (top) of characteristic wavelength 1 cm vs Ion Beam Sputtering 
(IBS) ripples (bottom) on Silicon. Water ripples were obtained by shaking a photography 
bucket filled with water (coutesy of R. Vida, Universidad Pontificia Comillas). Silicon 
ripples were obtained by bombarding Silicon with Ar + . The width of the observation 
window is 1 /im 2 (courtesy of L. Vazquez, Instituto de Ciencia de Materiales de Madrid- 
CSIC). 



ably, the outcome of these works invalidates the BH paradigm, since various 
morphological transitions are observed between unstructured and patterned 
surfaces as a function of both incidence angle 9 and ion energy E, confirming 
earlier observations at higher energies 17| and contradicting one of the main 
predictions of BH theory. Thus, ironically, while for the case of multiple- 
component systems differential sputtering and species segregation may have 
been recently identified 18M2l| as the main physical mechanisms behind the 
morphological instability, the a priori simpler case of single-component sys- 
tems still remains to be understood. 

In addition, none of the current phenomenological models that elaborate 
on the BH paradigm 22M24| are able to capture the complex morphological 
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features elucidated by the recent experiments on Si. In these theories, the 
dynamics of the target surface is described through an effective evolution 
equation in which the various physical mechanisms are combined ad-hoc, and 
defect dynamics are neglected. The same limitations when confronted with 
experiments occur even for improved physical models in which the evolution 
of the surface height is explicitly coupled to that of the density of defects 
whose transport is confined to a thin surface layer, akin to pattern formation 
on e.g. aeolian sand dunes, see e.g. 25l-l28| . and |7| for a recent overview. 



In order to overcome the limitations of previous continuum approaches to 
surface dynamics by IBS, we propose a generic framework based on general 
conservation laws (such as those governing mass or momentum), that allows 
to study systematically the role of the different mechanisms involved in the 
process. This approach, reminiscent of classical hydrodynamics, will allow 
us to identify the physics underlying the patterning of the eroded target 
surface. In summary, the main contribution of this work is to present a 
general consistent physical framework to study this system, that inherits the 
benefits of traditional fluid mechanics problems for which it has been proved 
highly successful in the past. 

2. Governing equations 

In the energy range considered, an impinging ion interacts mostly with 
the nuclear structure of the target, generating vacancies and interstitials. 
Ions and defects can diffuse at different rates, annihilate by recombination 
or form clusters. As borne out from Molecular Dynamics (MP) simulations, 
the overall effect of these processes is threefold: amorphizing [29] , stressing 



the target [30|, and displacing the mean atom position inside the material 



3 lj ] . Moreover, as also seen in experiments |3|], an amorphous layer forms 



rapidly [30|, |32j , its thickness and other mechanical properties like density, p, 
becoming stationary after a fluence of order 10 14 ions cm -2 (corresponding to 
a few seconds of irradiation for a typical ion flux of 10 13 to 10 15 ions cm -2 s _1 ). 
Such a small ion flux drives dynamics slowly, inducing (as in macroscopic solid 
flow for e.g. a glacier) a time scale separation between atomistic relaxation 
rates (~ 1 ps _1 for collision cascades or adatom hopping), and the ~ 10 s 
times in which significant morphological changes occur. This fact legitimates 
a continuum description. 
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2.1. Conservation of mass 

We are interested in the stages of the bombarding process in which the 
amorphous layer has already been formed. Thus, although its density value 
changes with respect to that for the pristine target, it does not evolve further 
after the initial transient just described. Mathematically, this fact can be 
expressed in the form 

V-V = 0, (1) 

where V is the velocity field of the fluidized layer, whose components are 
given by 

V = u\ + vj + iok. (2) 

In the language of fluid mechanics, the amorphous layer is said to be incom- 
pressible in the steady state. 

2.2. Conservation of momentum 



Conservation of momentum can be universally expressed as [33 



pDY/Dt = p(dtV + V • VV) = V ■ T, (3) 

where T is the stress tensor. Indeed, the physics of the IBS problem enters 
the constitutive equation for this tensor. 

Although a complete characterization of the stress tensor in the bom- 
barded material is still lacking, there are strong evidences that the amorphous 
layer can be considered as a highly viscous fluid, hence T can be written as 

Tij = -Pdij + n (d iUj + djUi) + T%, (4) 

where T?- contains all the terms that do not come from hydrostatic pressure, 
P, or viscous flow, and 1*1,2,3 = u,v,w are the components of the velocity 
field P). 

Moreover, because the radiation induced viscosity, /x, is around 15 orders 
of magnitude larger than the viscosity of water [9 ! ], we can safely drop the 
left hand side of Eq. ([3]). Thus, combining the conservation of mass with this 
assumption for the amorphous layer, we can write Eq. ([3]) simply as 

= V • T s + /iV 2 V - VP, (5) 

where we define b = V • T s as a body force acting in the bulk of the fluid 
layer. Actually, this force contains the relevant physics besides viscous flow. 
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It is effective in the sense that detailed knowledge of the stress induced be- 
neath the surface is not yet known experimentally. However, it is reasonable 
to expect that this force can be split into an amplitude, and an angular 
contribution, which is a function of the local angle of incidence. Mathe- 
matically, 



where 7 is the local slope of the surface and 6 is the incidence angle of 
the ion beam. Here, contains the coarse-grained information about the 
effect of the residual stress created in the target, due to ion-induced mass 
redistribution, and has dimensions of a gradient of stress. 

In three dimensions, and using standard spherical coordinates, we can 
write 



where <fi is the azimuthal angle (the angle in the plane of the target with 
respect to the plane of incidence of the ions) . Possible anisotropies can enter 
the dynamics from the difference in the values of the amplitudes Je and f' E . 
We are assuming that the collision cascade that produces the damage and, 
ultimately, the strain/stress that drives the amorphous layer, is oriented in 
the direction of the incoming ion (see Fig. 2). The fact that is the same 
in the equations for x and z comes from the fact that there are, actually, 
just two important directions: parallel and perpendicular to the ion flux. 
When we project the parallel contribution in those directions (x and z), the 
pre-f actor remains the same. 

This body force can be understood as an effective gravity term. Thus, in 
Fig. [2] we provide a cartoon to illustrate the analogy between ion induced solid 
flow and a genuine fluid flowing down an inclined plane. The two systems 
are mathematically equivalent once a rotation of angle 9 is performed in 
the laboratory reference frame. The main difference between both systems 
comes from the corresponding gravity fields. While for the standard fluid 
problem gravity is a constant, for the erosion system the effective gravity 
changes from point to point on the surface, its value depending strongly on 
the local surface orientation with respect to the ion beam. Thus, unlike the 
fluid problem, in the erosive case one can obtain a morphological instability 
for appropriate angles of incidence. 




(6) 




(7) 
(8) 
(9) 
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Figure 2: Analogy between IBS and fluid flow down a inclined plane. By rotating the 
frame of reference so that the ion impacts vertically, the IBS system can be seen as fluid 
flow down an inclined plane. The body force b is equivalent to a non-homogeneous gravity 
field that depends on the local slope value at the free interface. 
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2.3. Boundary conditions 

Equations ([I])-© need to be supplemented with proper boundary condi- 
tions involving the erosive terms and the stress at the interfaces. In the most 
general case, we have to specify the evolution for two important boundaries 
of the system: the amorphous- vacuum interface (the free surface), h^ a \ and 
the amorphous-crystalline interface, h^ c '. 

At the free surface, the balance of stress takes the form 1 3f 



n-T-n = —an + T^ xt \ (10) 
n-T-t = Va-t + T t (ext \ (11) 
n-T-q = Va-q + T^, (12) 

where n is the unit outward normal vector, t and q are two mutually per- 
pendicular vectors in the tangent plane to the surface, k is the mean surface 
curvature, and a is surface tension. Finally, T^ ext ^ are the projections along 
the corresponding directions of an external stress that is applied over the 
surface. 

Here we will consider that the surface tension is constant. However, 
for multicomponent materials where one species segregates at the surface, 
one can assume that there is a surface tension gradient that may affect the 



dynamics of the surface (the so-called Marangoni effect [33|, |34() 



Finally, a kinematic condition [33| at each interface leads to evolution 
equations for both hS a > and h^ c ' as 



Dt 

and 



W = jam, (13) 



—pf ~ W = Jer, (14) 

where, as defined above, w is the vertical component of the velocity field. 
The terms j am and j er stand for the rates of amorphization and erosion, 
respectively. Thus, Equations ( IT3|) and ( !T4|) express mathematically that the 
difference between the vertical velocity field and the actual motion of the 
interface is equal to the rate at which material is removed from both phases. 

For instance, the erosive (BH-type) mechanisms enter directly through 
j er , so that their contribution to the full dispersion relation will be additive, 
see section [31 For simplicity, we assume that the rate of amorphization is 
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equal to the rate of erosion, since in experiments both interfaces evolve at 
the same rates in the steady state. Moreover, for mathematical tractability, 
we assume that the crystalline- amorphous interface is flat. Under these as- 
sumptions, a final boundary condition is required, by which we specify that, 
at h( c \ the fluid moves with the same velocity as the crystalline phase; this 
is often referred to as the no-slip boundary condition |33j. 




3. Linear analysis 

The standard procedure to study the linear stability of a surface is devel- 
oped in three steps: 

1. Obtaining the flat (zeroth order) interface solution. 

2. Making a periodic infinitesimal perturbation of the solution found in 
the previous step, with wavelength A = 27r/g (q being the wavenumber). 

3. Calculating the amplification rate, u q , for every wavenumber q (also 
known as linear dispersion relation). 

For the sake of simplicity, we will detail these steps for the two dimensional 
case (namely, the interface is a line and not a proper surface). At the end 
of our derivation, we will quote the dispersion relation for the full three 
dimensional case. 

3.1. Flat interface solution 

We fix our reference frame on the free surface, so that the crystalline- 
amorphous interface is located at z = —d, where we denote by d the average 
thickness of the amorphous layer and by h(x) the fluctuations of the free 
surface. The flat interface condition is given simply by h(x) — 0. Introducing 
it in the equations for the velocity field, we find that the flat solution is further 
characterized by 




-f E zcoaO*(e), 
f E {d 2 -z 2 ) sin 6^(6) 



(15) 



ito(z) 



2/i 

0. 



(16) 
(17) 



w (z) 
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3.2. Linear perturbation 

To study the morphological stability of the flat interface, we look for 
solutions of the governing equations of the form 

h(x) = ee w «*- <9X , (18) 

for an arbitrary wavenumber q = 2tt/\. Similarly, we assume that the pres- 
sure and the velocity fields can be written as 

P(x,z) = P (z) +eP 1 (z)e u) " t - iqx , (19) 
u(x,z) = u (z) + eu l (z)e Ul > t - iqx , (20) 
w{x,z) = w (z) + ew l (z)e ul '' t - iqx , (21) 

and solve the governing equations perturbatively to linear order in the small 
parameter e. The precise form of the ensuing corrections to the flat interface 
solution, Pi(z), U\{z) and W\(z), is given by Eqs. (1A.ip - fjA.3p in Appendix A 
respectively. 

3.3. Dispersion relation 

Finally, the dispersion relation can be extracted from the kinematic con- 
dition. Specifically, from Eqs. (JHJ) and ( 1T8|) we get 

d t h = —u(x, 0)d x h + w(x, 0) + j er 
-t u q = iqu (0) + wi(0) + Q q , (22) 

where Cj q comes directly from the erosive contribution j er ; here is where, for 



instance, the classical theory of Bradley and Harper can enter [10 



Using Eqs. (|T5]) - (l2~T]) and denoting real part with a single prime and imag- 
inary part with a double prime, we find: 



q 2 a + f E d e (sm(9)V(6)) (-2dq + smh(2dq)) 

"J = 7 " x + "J (23) 

2qfi[l + 2d 2 q 2 + cosh(2dq) 



and 

[f E (d 2 q 2 (3 + 2dY + 



2q^[l + 2d 2 q 2 cosh(2dg) 
cosh(2dg)) sin(0)*(0) - 2 cos(#) (l + d 2 q 2 - 2 cosh(dg) + cosh(2dg) - 
2dq sinh(dg)) V'(9)\ ) + wj. (24) 
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3.3.1. Real part: Stability of the flat interface 

The real part of u q controls the stability of the flat interface: if u' q is pos- 
itive/negative, a periodic perturbation to the flat solution with wavenumber 
q will grow/decay exponentially as given by Eq. f fl8|) . Equation fl23|) gener- 
alizes important results from fluid dynamics, and has already been reported 
in the context of the IBS problem for more specific choices of the function \l/ 
0, [35|| . For instance, for /# — > 0, it reduces to Orchard's classic result 36 



on the flow of a viscous layer of arbitrary thickness on top of an immobile 



substrate. In particular, gcf > 1 then leads to Mullins' [37| celebrated rate of 



flattening when relaxation occurs through bulk viscous flow, u q = — crq/2fi, 



already invoked for some IBS experiments [38|, [39 . 

Generally, equation (123 p is non-polynomial, implying non-local effects in 
the dynamics of the amorphous layer. This non-locality may be important 
when other non-local effects are present in the dynamics (as redeposition, 
shadowing, or secondary ion scattering). Moreover, in our case the sign 
of u' may depend on the incidence angle, 6, thus explaining morphological 



transitions from flat to patterned structures at different values of 9 |13Hl6 



3.3.2. Imaginary part: in-plane propagation of the pattern 

In those cases in which an unstable wavelength is selected (namely, u' q > 
for some q), the imaginary part of u q provides information about the veloc- 
ity of lateral in-plane propagation of the pattern. Specifically, if the most 
unstable wavenumber is q c (namely, the one with the largest positive value 
of u' q ), then the pattern propagates with velocity 



V, 



dco" 



ripple 



dq 



<?=<?c 



(25) 



Again, the sign of V ripp i e may depend on the incidence angle; thus, depending 
on the choice of ^{0), the ripples can move in the same direction as or 
opposite to the ion beam. 

3.3.3. 'Shallow water" approximation (d/X <C 1 ) 

For large enough energies, and because the typical wavelengths of the 
patterns reported in the literature for energies around 1 keV are in the A w 
20 — 100 nm range, thus an order of magnitude larger than the amorphous 
layer thickness (d ~ 1—5 nm), we can assume that qd <C 1 and Taylor expand 



u' q . In fluid mechanics this is known as the shallow-water approximation [33 
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Thus, we get 



= 3^ q -^ q+ U <> (26) 



where it is understood that a similar expansion has been done in the erosive 
contribution uf . Note that, if we consider solid flow as the only pattern 
forming mechanism, we can simply drop this term. Since we are interested 
in probing the morphological consequences of viscous flow, we shall do this 
henceforth unless otherwise stated. 

Equation (125]) predicts that, if the function dg(^f(9) sm.9) is positive, the 
interface is stable (flat or rough, depending on the intensity of the thermal 
and beam flux fluctuations); on the contrary, if it is negative, a pattern 
appears with an angle dependent wavelength. 

At this point, note that our hydrodynamical formulation is a coarse 
grained description of the microscopic phenomena. Thus, it has been re- 
cently shown that the ion impacts produce a response in the target referred 



to as a prompt regime 40] . In this reference, starting from atomistic molec- 
ular dynamics (MD) simulations, the overall effect of the ion impacts on the 
evolution of the surface height is summarized into the moments, M^, of a 
so-called crater function, that eventually contribute to a continuum evolu- 
tion equation for the surface height. It is worth noting that we can map the 



approach 40j with our present formalism through 

#(0) = M£Vtau(0). (27) 

In this respect, we could incorporate MD information into our hydrodynam- 
ical description through our function Let us note incidentally that, as 



suggested by the simulations in |40|, |4l| , purely erosive contributions to the 



evolution of the surface height are numerically negligible, justifying our ne- 
glect of ZJ' q above. 

Patterns generated in the unstable parameter region in which Eq. ( )26l) 
has a positive maximum for q = q c , have a characteristic wavelength given 
by 

2tt / 2a N 1/2 



Xc -^- 2n {-f E d e (mM0))) ■ (28) 

For instance, if one assumes that the angular correction is ^(6) = cos 6 
(namely, that the body force is proportional to the local flux at the surface), 
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then Eq. f[2"6"j) predicts a transition from flat to patterned surfaces at 9 = 45° 
for all ion energies, as has been observed in experiments of Ar + on Si at 
energies between 250 eV and 1 keV 14|. Similarly, for this choice of the 
velocity of ripple is given by 

Vrippie = sm 29 II H — I . (29) 

To be more quantitative, we need to determine the values of the main 
parameters of the theory. As mentioned above, the new parameter Je is 
related to the stress gradient created across the amorphous layer, induced 
by the collisional processes. There is a large experimental uncertainty in the 
measurement of stress, a property that depends moreover on the energy of the 
incident ions. Different authors report different values for ion-induced stress 
on Si, in the range 200 MPa to 1.62 GPa M,M. We are interested in the 
order of magnitude, hence we take the geometrical mean (which accounts for 
the average order of magnitude) of both extreme values, thus, approximately 
569 MPa. This corresponds to = 0.424 kg nm _2 s~ 2 for a layer of d = 3 
nm depth. In Fig. 3, this value corresponds to the blue solid line, while the 
limiting values of 200 MPa and 1.62 GPa are used to plot the red dashed 
lines, which we take as our confidence interval. More accurate measurements 
of stress would provide improved estimations. Additional parameters have 
not been measured experimentally yet. However, we can provide indirect 
estimates for their values when these are not explicitly available, as shown 
in Table 1. With these parameters at hand, Eq. (128]) and experimental data 
are seen to agree quite closely in Fig. [3j 



Parameter 


Value 


Reference 


Surface tension (a) 


1.34 J/m 2 


[42] 


Amorphous layer thickness (250 eV) 


3 nm 


[13, 14] 


Viscosity (fj,) 


6 x 10 8 Pa s 


[40]* 


Stress range 


200 MPa-1.62 GPa 


[14,30] 



Table 1: Numerical values of the parameters used in the main text. *Mayr et al. [9( report 
the value n = 10 9 Pa s, which also leads to good agreement with the results presented in 
this work. 



The dependence of the pattern properties on other parameters such as 
energy, flux, etc., that enter the amplitude will be the subject of future 
work. 
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Figure 3: X c (0) at the type II transition. Experimental data from 13j, [l4[ are shown as 
black symbols, the green vertical line is the transition angle, and the blue solid line is the 
present theoretical prediction from Eq. (|28|) : the red dashed lines provide our confidence 
interval, due to uncertainties in parameter values corresponding to experimental data. A 
is nm and 9 in degrees. 



3.4- Dispersion relation for the three dimensional case 

For completeness, we include here the dispersion relation for the three 
dimensional case. Now, the flat interface is a plane, and deviations to linear 
order take the form 

h(x,y) = ee^- iqoc , (30) 

with q = q x i+q y i and x = xi+yj. In the "shallow-water" limit, the dispersion 
relation is now given by 



f' E cos 9V(6)d 3 2 f E d e (sm9^(9))d 3 2 ad 3 2 2 2 , 



for a general angular dependence \I> of the body force and general (anisotropic) 
forms for the strain induced by the ion through the amplitudes f E and f' E . 



4. Connection with other continuum frameworks 

4-1- Two-field "hydrodynamic" theory 

Despite the fact that we have not considered the dynamics of the amorphous- 
crystalline interface, here we want to note that, in the shallow-water limit, 
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one can reduce the present full hydro dynamical formulation, making contact 
with previous theories of erosion, specifically with the two-field model ap- 

see an overview in 



roach pioneered in |43| and further developed in |25H28 



Within this approach (see Ref. ji?} for details), the dynamics of the den- 
sity R of moving species and that of the free surface height, h, are explicitly 
described, within the implicit assumption that the thickness of the surface 
layer where material transport takes place is negligible. The evolution of 
these two fields can be expressed through proper rate equations, borrowed 
from the field of pattern formation on aeolian sand dunes. Specifically, 

d t R = (l-(f>)T ex -T ad + DV 2 R-V-VR, (32) 
d t h = -T ex + T ad , (33) 

where (1—0) measures the fraction of atoms that, having been dislodged from 
their equilibrium positions, are actually sputtered. Here, T ex and T ad are, 
respectively, rates of atom excavation from and addition to the (crystalline) 
immobile bulk. 

In the hydrodynamical framework proposed in the present work, the kine- 
matic conditions ([TBI and f[T4"|) can be combined and reinterpreted within the 
two-field model approach. Thus, if we define 

R = h^-h^ c \ (34) 

then 

DR 

— = d t R + V-VR = j er - j am . (35) 
Finally, we just need to define 

j er = -0T ex + DV 2 R, (36) 

jam Tex Tad- (37) 

to recover Eq. (132]) exactly. 

4-2. Lubrication theory 

The hydrodynamical nature of the model considered in Sections [21 [3] sug- 
gests a direct comparison with the so-called lubrication theory of fluids. This 
theory is valid in the shallow- water approximation which is, in fact, the one 
relevant to IBS experiments. 
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Following Oron et al. 33(, we define a small parameter e = qd. Addition- 
ally, we assume that not only the slopes are small but also the time-scales 
related to the motion of the fluid are slow (with respect to the motion of 
the very same surface). Mathematically, the following change of variables is 
defined: 

X = ex/d, Z = z/d, U = u/U , W = w/eU and T = eU t/d. (38) 

In order to capture the essential ingredients that affect the motion of the fluid, 
we rescale the stress tensor and the corresponding values at the boundary 
Thus, 

d ^ d d 

7~xx j j xx j T X z j j xz i 7~zz j j zz 5 v^^J 

eU U EUq 

and 

p = £ j^p, ^ { : xt) = £ 4 T ^ ext) . ^ xt) = tj T t ext) ■ ( 4 °) 

U Uq Uo 

Finally, the surface tension in the proper units is redefined as 33( a = e 3 a. 



Introducing this change of variables into the governing equations and bound- 
ary conditions, and keeping the dominant contribution to order 0(e°), we 
find 

d x U + d z W = 0, (41) 
-d x P + dzZ xz = B x , (42) 
-d z P = -B z , (43) 
P\z=h(x,t) = -ad 2 x H-^\ (44) 
^xz\z=h(x,t) = T,[ ext) +d x a. (45) 
Integrating this system, one arrives to 

r H(X,T) 

d T H + d x / U(X,Z,T)dZ = -J er , (46) 

Jo 

where J er is the contribution from erosion in the scaled variables. 

In order to solve ( )42l) -f j45l) . one needs some closure relation between the 
stress and the velocity. In our problem, we have split the stress into a body 
force and a (Newtonian) viscous tensor. For such a case, one arrives to a 
closed non-linear equation of the form 

d T H = -d x \^(-B x + B z Hx + aHxxx + dx^ n ext) ) 
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HUB Z rf 7 (~{e*t)^ da \ H 2 



J P: 



with 



B x 

P(X,Z,T) 
B z 



- 7) sin 9, 
-B Z {Z -H) + ad 2 x H - T^ xt \ 



Uo 
d 2 



f E ^{o-i) cose. 



(47) 
(48) 
(49) 



Hence, with iff (9) = cos9 and F E = f E U /d 2 , 

~F R H 3 f-cos9 sin 9 - cos 29 H x + sin 9 cos 9H 



d r H 



-d 



x 



X 



H — ^-crHxxx 

J pr • 



(i + H 2 x y/ 2 

F E H XX H 4 sin 9 cos 9- cos 2 9H X 



[1 + H 2 X )V 2 



(50) 
(51) 



Within our framework, an alternative equation can also be derived that fea- 
tures the same real part of the linear dispersion relation, setting B X z = 0, 
Y*(ext) _ ^ j- n Y>n Xt ^ = n ■ T ■ n, with T as derived from Ref. for 
a two-dimensional coordinate system aligned with the ion beam, 

— cos 9 2 + v sin 2 9 — (1 + v) cos # sin £ 

— (1 + v) cos # sin # v cos 2 9 — sin 2 9 



r 



F' E 



(52) 



where v 
d T H = 



1 for a two-dimensional incompressible fluid. Thus 

-d x 
H 3 



F' E H 2 



{H% - 1) sin 29 + 2 cos 29H X 



1 + H 2 



x 



+— uH xxx + F'd x 



1 - H x )cos29 + 2sin29H 



x 



1 + H* 

with Fe 7^ F' E < for compressive stress (the experimental case). 

Linearizing Eq. ( 14 7p and neglecting the effects of external stresses and the 
Marangoni term (d x a), we reproduce the shallow water dispersion relation 
given by Eq. (126]) . The interest of this lubrication approach is that, provided 
an appropriate closure relation between the stress tensor and the velocity field 
is put into Eq. f|42|) . we obtain a non-linear equation for the surface height 
through Eq. P6l) . given Eq. P2|) . Thus, nonlinear effects can be assessed 
that go beyond the analysis performed in Section [3J 
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5. Conclusions 



In this work we have introduced a physical framework for surface pattern 
formation by IBS that is inspired in classic Fluid Mechanics. The strength 
of the theory is that it is based on general conservation laws, and that all the 
physics of the erosive process can be included into appropriate constitutive 
equations for the stress (body force) and boundary conditions. This approach 
allows to properly account for different mechanisms, such as viscous flow, 
erosion, or ion-induced strain (damage). 

We have performed a linear stability analysis of the governing equations 
and shown that the forces related to the stress produced by the ion impacts 
indeed influence the morphological stability of the interface. For a specific 
choice of the angular dependence of the body force, \1>, we have seen that 
some of the complexities of the morphological diagram for Si that have been 
unveiled recently, can be accounted for by our description. In this connection, 
purely erosive contributions indeed seem to play a minor role in the surface 
dynamics as compared with viscous flow. 

Finally, the proposed hydrodynamic framework can be seen to generalize 
previous two-field formulations of IBS, to the case of a flowing layer of ar- 
bitrary thickness. Such phenomenological models have proved successful to 
obtain an effective equation describing the non-linear dynamics of the surface, 



reaching in some cases quantitative [4J and semiquantitative 45j agreement 



with experiments. On the other hand, our present framework does allow for a 
systematic study of nonlinear effects through a lubrication theory approach. 

The specific details of the constitutive equations that are appropriate for 
different experimental conditions will be the aim of future work. Also the 



general non-linear theory will be addressed through standard procedures [33 
that should be able to account for nonlinear phenomena such as pattern 
coarsening, saturation of the roughness, or defect creation and/or annihila- 
tion. 
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Appendix A. Solutions of the linearized governing equations 

Expressions for Pi (z), U\(z) and W\(z) in Eqs. (|T9j) to (I2T]) of the main 
text are 



1 



-q(2d+z) 



2\l + 2d 2 q 2 + cosh(2dq) 

q 2 (l + e 2dq (l + 2dq + e 2qz (l + e 2dq - 2dq^j \ \ a + 2e q(2d+z) (f E cos(9) (cosh(qz) + 
cosh(g(2<i + z)) — 2dgsinh(gz)) + i/g sm(0)(2G?g cosh(g;z) + sinh(gz) + 
sinh(g(2ci + (0) + (i/e (l - e d " (2 + 2e Mg - 2e 2qz + e 3dq+2qz - 

Adq + e dq {-l + 2dq) - 2e 2q{d+z \l + 2dq)e q{d+2z) (1 + 2dg))) cos(#) - 

/ E (l - e qz + e 2q{2d+z) - e q{4d+z) + e 2qi - d+z \l - 2dq) + e 2dq (l + 2dq) - 

2e q(2d+z) ^ + 2d 2 q 2 ^ sin(0))^'(0) , (A.l) 



Ui(« 



"?2 



2qfi (l + e 4d ? + e 2d i (2 + Ad?q 2 
3 ^ _ e 2,(2d +2 ) 2 + e 2,(d+,) + 2dg(d + *)) + e M<? (^ + 2dg(d + z))) a + 

f E (iq(z - e 2q{2d+z) z + e 2q{d+z \-z + 2dq{d + z)) + e 2 ^ + 2dq(d + z))) co 

( - 1 + qz + e 2qi2d+z \l + qz) + e 2d<z (-l + q(z - 2d(-l + q{d + *)))) + 

e 29(d+2) (l + q(z + 2d(l + q(d + 2))))) sin(0)W(0) + 

f E f(l-qz + 2e 3dq q{d + z) + 2e q(d+2z) q{d + z) - e 2qi - 2d+z \l + qz) 
-2e dq q{-d -z + 2dqz) + 2e 3dq+2qz q(d + z + 2dqz) + e 2dq (l + 2dq{-\ + dq) 
+q(-l + 2dq)z) - e 2q(d+z) {l + q(z + 2d(l + q{d + *))))) cos(#) 



-iq{z 



^2q(2d+z) 



z + e 2q{d+z) (-z + 2dq(d + z)) + e 2dq (z + 



2dq(d + z)))sm(6))V(6) 



(A.2) 
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and 



wiiz) = ; ; ^e~ qz 



2g/i (l + e Mc > + e M <? (2 + 4d 2 g 2 ) ) 
q 2 (l + qz + e 2q{2d+z \-l + qz) + e 2q{d+z \-l + q(z - 2d(-l + q(d + z)))) + 
e 2d<z (l + g(z + 2d(l + q(d + z)))))<7 + + + e 2 ^ M+2 )(-l + gz) + 

e 2 q (d+z) + ^ _ 2 ^_ 1 + g ( d + + + ^ + 2d{l + ^ + z))))) cos(#) 

+if E q(~ z + e 2q(2d+z) z + e 2dq (-z + 2dq(d + z)) + e 2g(d+2) (,2 + 

2dg(d + z)fj sin(0))*(0) + (if E (2e qz + 2e <?(4d+2) + e ?(2d+2) (4 + 8rf 2 g 2 ) + 

qz _ e 2q{2d+z) qz + 2e d9 (-l + + g(-l + 2dq)z) + 2e 9(d+22) (-l + q(d + z)) - 
2e Mg (l + q(d + z)) - e 2dq q{-z + 2dq[d + z)) - e 2q{d+z ^q(z + 2dq[d + z)) + 
2e 3d 9 +2,z(_ 1 + g2 + dq (_ x + 2gf ^^ cog ^) + / Bsin (0)( _ i _ qz + eMM+z)^ _ qz) + 

2e<? (2d+ 2 ) ^ _ g ^ 2rf + ^ cosh ( g2 ) + (l + 2rfg 2 (rf + *)) sinh(gz)) ) ) #'(0) . (A.3) 
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